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INTRODUCTION 


The  piecewise  linear  analysis  option  of  the  NASTRAN  program  provides  an 
algorithm  for  solving  nonlinear  problems  in  material  plasticity. 1  However, 
the  usefulness  of  this  option  is  quite  limited  because  only  a  few  elements 
have  been  implemented.  These  include  rod,  tube,  and  bar  elements  for  one¬ 
dimensional  problems  and  plate  elements  for  two-dimensional  plane  stress  prob¬ 
lems.  In  a  recent  paper, 2  the  implementation  of  a  trapezoidal  ring  element  in 
NASTRAN  for  elastic-plastic  analysis  was  described  and  two  test  problems  of 
rotational  symmetry  were  solved.  The  first  is  an  infinitely  long  tube  under 
uniform  internal  pressure.  The  NASTRAN  results  are  in  excellent  agreement 
with  an  exact  solution  based  on  the  finite-difference  approach.^  The  second 
problem  is  a  thick-walled  cylinder  of  finite  length  loaded  over  part  of  its 
inner  surface.  The  NASTRAN  results  are  in  good  agreement  with  another  finite 
element  solution.^ 

In  the  present  report,  two  elastic-plastic  problems  of  rotational  symmetry 
are  solved  using  triangular  ring  elements  for  the  finite  element  models.  The 
implementation  of  a  triangular  ring  element  in  NASTRAN  for  elastic-plastic 
analysis  follows  the  same  procedures  in  reference  2  for  a  trapezoidal  ring 

1mNASTRAN  Theoretical  Manual,”  NASA  SP-221(01),  1972. 

^Chen,  P.  C.  T.  and  O'Hara,  G.  P. ,  "Implementation  of  a  Trapezoidal  Ring 
Element  in  NASTRAN  for  Elastic-Plastic  Analysis,"  NASA  CP-2131,  Proceedings 
of  the  Eighth  NASTRAN  User’s  Colloquium,  pp.  101-112,  October  1979. 

3Chen,  P.  C.  T.,  "A  Finite-Difference  Approach  to  Axisymmetric  Plane-Strain 
Problems  Beyond  the  Elastic  Limit, “  Transactions  of  the  Twenty-Fifth 
Conference  of  Army  Mathematicians,  pp.  455-465,  January  1980. 

^Chen,  P.  C.  T.,  "Numerical  Solution  of  Gun  Tube  Problems  in  the  Elastic- 
Plastic  Range,"  Proceedings  of  the  1977  Army  Numerical  Analysis  and  Computer 
Conference,  pp.  423-439,  November  1977. 
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element.  In  order  to  test  the  accuracy  of  the  implementation,  the  plane- 
strain  problem  of  a  part ially-plast ic  thick-walled  cylinder  under  internal 
pressure  is  solved  again  and  compared  with  the  finite-difference  solution.^ 

The  second  demonstrative  problem  in  using  NASTRAN  triangular  ring  elements  is 
an  overloaded  thread  problem  for  the  British  Standard  Buttress.  A  detailed 
discussion  of  this  problem  in  the  elastic  range  of  loading  was  reported  in 
reference  5.  In  the  present  report,  uniform  pressure  distribution  on  the  pri¬ 
mary  bearing  surface  of  the  thread  is  assumed  and  the  load  is  applied  in 
increments  beyond  the  elastic  limit.  The  maximum  axial  and  principal  stresses 
are  located  and  their  values  are  determined  as  functions  of  loadings. 

TRIANGULAR  RING  ELEMENTS 

The  incremental  displacement  field  employed  for  the  triangular  ring 
element  is 

Au(r , z)  *=  +  $2r  +  $3Z  >  (1) 

Aw(r  ,z)  =  +  35r  +  06Z  .  (2) 

The  transformation  from  grid  point  coordinates  to  generalized  coordinates  is 

{e}  =  [r0q]  {Aq}  (3) 

where 

{ q} T  =  Awl  ,  Au2,  Aw2  ,  *  Au3  ,  Aw3J  ,  (4) 

{3}t  =  | _ 3 1 ,  &2>  33,  04,  05,  36^,  (5) 


^Chen,  P.  C.  T. ,  ”A  Finite-Difference  Approach  to  Axisymmetric  Plane-Strain 
Problems  Beyond  the  Elastic  Limit,'1  Transactions  of  the  Twenty-Fifth 
Conference  of  Army  Mathematicians,  pp.  455-465,  January  1980. 

^O’Hara,  G.  P.,  "Stress  Concentrations  in  Screw  Threads,"  NASA  CP-2131, 
Proceedings  of  the  Eighth  NASTRAN  User’s  Colloquium,  pp.  65-77,  October  1979. 
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and  the  elements  of  the  inverse  of  the  transformation  matrix  [rgq^1  are  the 
coefficients  of  the  3 1  s  in  equations  (1)  and  (2),  evaluated  at  the  corners  of 
the  element. 

The  stiffness  matrix  is  formed  in  the  same  manner  as  that  for  the 
anisotropic  elastic  element.  The  final  form  referred  to  grid  coordinates  is 

[k]  =  [rgq]T  [k]  [rBq]  ,  (6) 

where 

[k]  =  2r/r [ b] t  [d]  [b]  dzdr  .  (7) 

The  [ b]  matrix  is  the  same  as  the  elastic  case,  but  now  it  expresses  the 
incremental  strains  in  terms  of  generalized  coordinates 

{At}  =  [b]  {s}  .  (8) 

The  [d]  matrix  which  relates  the  incremental  stresses  to  the  incremental 
strains,  i.e., 

{Ao}  =  [d]  {Ae}  ,  (9) 

is  the  same  as  that  for  a  trapezoidal  ring  element  presented  in  reference  2. 

In  developing  the  NASTRAN  program  for  strain-hardening  materials,  we  calculate 
[  D] -1  and  obtain  its  inverse  [ d]  numerically.  For  ideally-plastic  materials, 
this  procedure  fails  and  we  should  calculate  [d]  directly  using  the  closed 
form. 


^Chen,  P.  C.T.  and  O'Hara,  G.  P.,  "Implementation  of  a  Trapezoidal  Ring 
Element  in  NASTRAN  for  Elastic-Plastic  Analysis,"  NASA  CP-2131,  Proceedings 
of  the  Eighth  NASTRAN  User's  Colloquium,  pp.  101-112,  October  1979. 
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NASTRAN  IMPLEMENTATION 


The  implementation  of  a  triangular  ring  element  in  NASTRAN  for  elastic- 
plastic  analysis  follows  the  same  procedures  for  a  trapezoidal  ring 
element.2 * * *,6  Changes  were  required  in  the  functional  modules  PLA1 ,  PLA3 ,  and 
PLA4 ,  which  included  the  writing  of  seven  new  subroutines.  The  changes  in 
PLA1  allow  this  module  to  identify  the  new  element  as  a  member  of  the  piece¬ 
wise  linear  element  set  and  properly  initialize  the  nonlinear  Element  Summary 
and  Element  Connection  Property  Tables.  Three  element  stress  recovery  sub¬ 
routines  were  added  to  PLA3 :  PSIARG,  a  driver  for  stress  data  recovery; 

PSRIR1  and  PSRIR2 ,  phase  I  and  II  stress  recovery  routines.  Element  stiffness 
calculations  in  PLA4  required  four  new  subroutines:  PKIARG,  a  driver  for  non¬ 
linear  triangular  ring  elements  in  PLA4 ;  PKRIR1  and  PKRIR2 ,  stress  recovery 
routines  which  generate  stresses  for  the  computation  of  the  nonlinear  material 
matrix;  and  PKRIRS,  the  stiffness  matrix  generation  routine  for  nonlinear  trL- 
angular  ring  elements. 

PROGRAM  EVALUATION 

In  order  to  evaluate  the  accuracy  of  the  computed  code,  the  plane-strain 
problem  of  an  elastic-plastic  thick-walled  cylinder  under  uniform  internal 
pressure  is  solved  again  and  compared  with  the  finite-difference  solution.2 

2Chen,  P.  C.  T.  and  O'Hara,  G.  P.,  "Implementation  of  a  Trapezoidal  Ring 
Element  in  NASTRAN  for  Elastic-Plastic  Analysis,"  NASA  CP-2131,  Proceedings 
of  the  Eighth  NASTRAN  User's  Colloquium,  pp.  101-112,  October  1979. 

Chen,  P.  C.  T.,  "A  Finite-Difference  Approach  to  Axisymmetric  Plane-Strain 

Problems  Beyond  the  Elastic  Limit,"  Transactions  of  the  Twenty-Fifth 

Conference  of  Army  Mathematicians,  pp.  455-465,  January  1980. 

6"NASTRAN  Programmer's  Manual,"  NASA  SP-223(01),  1972. 
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The  tube  of  outside  radius  2  inches  and  inside  radius  1  inch  is  divided  into 
10  equal  intervals  and  each  interval  consists  of  2  triangular  ring  elements. 
The  material  constants  are  E  =  30  x  106  psi,  v  =  0.3,  0O  =  1.5  x  105  psi  and 
the  effective  stress— strain  curve  is  represented  by  three  line  segments 
connecting  the  four  points  in  the  (e  ,a)  plane,  (e,o/oQ)  =  (0.0, 0.0), 
(0.005,1.0),  (0.055,1.5),  (0.1, 1.5).  Twenty-three  load  factors  are  chosen: 
P/°o  =  0.4323,  0.4738,  0.5125,  0.5484,  0.5818,  0.6128,  0.6415,  0.6681,  0.6925, 
0.7150,  0.7356,  0.7545,  0.7716,  0.7871,  0.8011,  0.8135,  0.8245,  0.8341, 

0.8423,  0.8493,  0.8550,  0.863,  0.87.  The  numerical  results  based  on  the 
NASTRAN  program  have  been  obtained.  For  this  problem,  an  exact  solution  based 
on  a  new  finite— difference  approach3  can  be  used  to  assess  the  accuracy  of  the 
NASTRAN  code.  Some  of  the  results  for  the  displacements  and  stresses  are 
Presented  graphically  in  Figures  1  and  2.  The  radial  displacements  at  the 
inside  as  well  as  outside  surface  are  shown  in  Figure  1  as  functions  of  inter¬ 
nal  pressure.  Figure  2  shows  the  distributions  of  radial,  tangential  and  axi— 
^1  stress  components  in  a  part ially— plast ic  tube  when  the  pressure  is  p  = 
0.7356  o0.  As  demonstrated  in  Figures  1  and  2,  the  NASTRAN  results  are  in 
very  good  agreement  with  those  based  on  the  finite-difference  approach.3 


3Chen,  P.  C.  T. ,  A  Finite-Difference  Approach  to  Axisymmetric  Plane-Strain 
Problems  Beyond  the  Elastic  Limit,"  Transactions  of  the  Twenty-Fifth 
Conference  of  Army  Mathematicians,  pp.  455-465,  January  1980. 
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THREAD  PROBLEM 


As  an  application  of  using  NASTRAN  triangular  ring  elements  for  two- 
dimensional  problems,  an  overloaded  thread  problem  for  the  British  Standard 
Buttress  is  examined.  A  finite  element  model  is  shown  in  Figure  3.  A  detail¬ 
ed  discussion  of  this  problem  in  the  elastic  range  of  loading  was  reported  ir 
reference  5.  In  the  present  report,  uniform  pressure  distribution  (P)  on  the 
primary  bearing  surface  of  the  thread  is  assumed  and  the  load  is  applied  in 
ten  unequal  increments  beyond  the  elastic  limit.  The  elastic  constants  are  E 
=  25000  Ksi,  v  =  0.25,  and  the  effective  stress-strain  curve  is  represented  by 
five  line  segments  connecting  the  six  points  in  the  (e-a)  plane,  (e-a  in  Ksi) 

=  (0.0,  0.0),  (0.0048,120),  (0.0122,180),  (0.0167,200),  (0.0918,210), 
(0.25,210).  The  sides  A-B  and  C-D  are  constrained  in  the  axial  direction  and 
the  side  B-C,  in  the  axial  and  radial  directions. 

The  elastic  problem  is  solved  first  and  the  upper  limit  of  the  elastic 
loading  (p*)  is  found  to  be  54.65  Ksi.  Eleven  load  factors  are  chosen  as  fol¬ 
lows:  p  =  54.65,  60.59,  66.22,  71.54,  76.54,  81.23,  85.61,  89.68,  93.43, 
96.87,  100.0  Ksi.  The  numerical  results  based  on  the  NASTRAN  code  have  been 
obtained  and  the  total  CPU  time  is  59  minutes  on  IBM  360  Model  44.  The  plas¬ 
tic  zone  at  the  maximum  load  is  indicated  by  the  shaded  area  in  Figure  3. 

Some  of  the  stress  results  are  shown  graphically  in  Figures  4  through  6.  Of 
particular  interest  is  the  region  along  side  D— E— F .  The  axial  stress  (az)  and 

^O'Hara,  G.  P.,  "Stress  Concentrations  in  Screw  Threads,"  NASA  CP-2131, 
Proceedings  of  the  Eighth  NASTRAN  User's  Colloquium,  pp.  65-77,  October  1979. 
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major  principal  stress  (a^)  in  the  boundary  elements  along  side  D-E-F  are 
shown  in  Figures  4  and  5,  for  the  first  and  last  load  factor,  respectively. 
The  maximum  fillet  stress  (a^)  occurs  near  E,  while  the  maximum  axial  stress 
always  occurs  at  D.  Finally,  the  maximum  fillet  stress  and  axial  stress  are 
plotted  as  functions  of  contact  pressure  in  Figure  6. 

CONCLUSION 

An  elastic-plastic  triangular  ring  element  has  been  implemented  in  the 
NASTRAN  computer  program.  Its  accuracy  has  been  evaluated  by  solving  a 
simpler  problem  for  which  an  exact  solution  is  available.  Its  application  to 
more  general  problems  has  also  been  demonstrated  by  solving  an  overloaded 
thread  problem. 
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Figure  1.  Radial  displacement  in  a  pressurized  tube. 
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Figure  2.  Stresses  in  a  pressurized  tube. 
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Figure  3.  Finite  element  model  for  the  thread  problem. 
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Figure  4.  Stresses  in  boundary  elements  (side  D-E-F)  at  p  =  54,65  Ksi. 


12 


*10  ksi 


Figure  5.  Stresses  in  boundary  elements  (side  D-E-F)  at  p  =  100  Ksi. 
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MAX.  AXIAL  &  FILLET  STRESS,  10  ksi 
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6.  Maximum  axial  and  fillet  stresses  as  functions  of  contact  pressure. 
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WARREN  MICHIGAN  4 '8090 


CCMMANDER 

US  MILITARY  ACADEMY 

ATTN:  CHMN,  MECH  ENGR  DEPT  1 

WEST  POINT,  NY  10996 

CCMMANDER 
REDSTONE  ARSENAL 

A  TIN:  DRSMI-RB  2 

-RRS  1 

-RSM  1 

ALABAMA  35809 


CCMMANDER 

ROCK  ISLAND  ARSENAL 

ATTN:  SARRI-ENM  (MAT  SCI  DIV)  1 

ROCK  ISLAND,  IL  61202 

CCMMANDER 

HQ,  US  ARMY  AVN  SCH 

ATTN:  OFC  OF  THE  LIBRARIAN  1 

FT  RUCKER,  ALABAMA  36362 

CCMMANDER 

US  ARMY  FGN  SCIENCE  &  TECH  CEN 
ATTN:  DRXST-SD  1 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 

CCMMANDER 

US  ARMY  MATERIALS  &  MECHANICS 
RESEARCH  CENTER 

ATTN:  TECH  LIB  -  DRXMR-PL  2 

WATERTOWN,  MASS  02172 


PLEASE  NOTIFY  CCMMANDER,  ARRADCCM,  ATTN:  BENET  WEAPONS  LABORATORY, 
DRDAR-LCB-TL,  WATERVLIET  ARSENAL,  WATEKVLIET,  N.Y.  12189,  OF  ANY 
REQUIRED  CHANGES. 
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COMMANDER 

US  ARMY  R]iS LARCH  OFFICE- 
P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709 
CORIANDER 

US  ARMY  HARFIY  DIAMOND  LAB 
ATTN:  TECH  LIB 
2800  POWDER  MILL  ROAD 
ADELPHIA,  MI  20783 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENG  ACT 

ATTN:  DRXPE-MT 

ROCK  ISLAND,  IL  61201 

CHIEF,  MATERIALS  BRANCH 
US  ARMY  R&S  GROUP,,  EUR 
BOX  65,  FPO  N.Y.  09510 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CEN 
ATTN:  CHIEF,  MAT  :SCIENCE  DIV 
DAHLGREN,  VA  22448 

DIRECTOR 

US  NAVAL  RESEARCH  LAB 
ATTN:  DIR,  MECH  DIV 

CODE  26-27  (DOC  LIB) 
WASHINGTON,  D.  C.  20375 

NASA  SCIENTIFIC  &  TECH  INFO  FAC. 

P.  0.  BOX  S757,  ATTN:  ACQ.  BR 
BA LT IM 0RE/WA3H INGTCN  INTL  AIRPORT 
MARYLAND  21 240 


COMMANDER 

DEFENSE  TECHNICAL  INFO  CENTER 
1  ATTN:  DTLA-TCA 

CAMERON  STATION 
ALEXANDRIA,  VA  22314 

META LS  &  CERAMICS  INFO  CEN 

1  BATTELLE  COLUMBUS  LAB 

505  KING  AVE 
COLUMBUS,  OHIO  43201 

MECHANICAL  PROPERTIES  DATA  CTR 
BATTELLE  COLUMBUS  LAB 

1  505  XING  AVE 

COLUMBUS,  OHIO  43201 

MATERIEL  SYSTEMS  ANALYSIS  ACTV 
1  ATTN:  DRXSY-MP 

ABERDEEN  PROVING  GROUND 
MARYLAND  21005 


1 

1 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARRADCCM,  ATTN:  BENET  WEAPONS  LABORATORY, 
DRTAF- LCB-TL,  WATERVLIET  ARSENAL,  WATERVLIST,  N.Y.  12189,  OF  ANY 
REQUIRED  CHANGES. 


NO.  OF 
COPIES 


12 


1 


1 


